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Abstract 

In this paper, we study the following knapsack problem: Given a list of squares with profits, we 
Q\ ' are requested to pack a sublist of them into a rectangular bin (not a unit square bin) to make profits 

in the bin as large as possible. We first observe there is a Polynomial Time Approximation Scheme 
(PTAS) for the problem of packing weighted squares into rectangular bins with large resources, 
^ , then apply the PTAS to the problem of packing squares with profits into a rectangular bin and get 

^—^ ' a | + e approximation algorithm. 

lHj' 1 Introduction 
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In this paper we study the two-dimensional generalization for Knapsack: we are given a set of squares, 
each of which is associated with a profit. The goal is to pack a subset of the squares (items) into a 
rectangle (bin) to maximize the total profit packed. The problem is NP-hard in the strong sense even 
^ \ if each item is an unweighted square (i.e., its profit is equal to its area) p2]. A little surprisingly, the 
research for approximation algorithm has started quite recently: Jansen and Zhang [12], Caprara and 
Monaci [2], Harren [IT] etc. 

Related Work There are many literatures on rectangle packing and square packing. For a 
two dimensional knapsack problem in which a subset of a given set of rectangles are packed into a 
given rectangular bin to maximize the total profits in the bin. Jansen and Zhang proposed 2 + e 
^ approximation algorithm [12]. When all items are squares and their profits are equal to their areas, 
Fishkin, Gerber, Jansen and Solis-Oba [8] presented a PTAS, which was also obtained by Han, Iwama 
and Zhang independently [ID] . Jansen and Zhang [13] proposed a PTAS for packing squares into 
a rectangular bin to maximize the number of squares packed in the bin [13] . Harren [IT] proposed 
| + e approximation algorithm for packing squares into a unit square bin. But his algorithm is not 
applicable to pack squares into a rectangular bin since his algorithm requires that every side of the 
bin must have the same length. Fishkin, Gerber and Jansen [7] obtained a (1 — e)-approximation 
algorithm for packing a set of rectangles with profits into a large resource bin with width 1 and height 
larger than (1/e) . 

Another related work is 2 dimensional bin packing problem in which all rectangles have to be 
packed into a unit square bin to minimize the number of bins required. When all items are squares, 
Ferreira et al. [6] gave an approximation algorithm with asymptotic worst-case ratio bounded above 
by 1.988. Kohayakawa et al. [16] and Seiden and van Stee |19j independently obtained approximation 
algorithms with asymptotic worst-case ratio of at most 14/9 + e (for any e > 0). These results were 
recently improved by Bansal, Correa, Kenyon and Sviridenko |lj. They proposed an asymptotic PTAS 
for packing d-dimensional cubes into the minimum number of unit cubes. For the online case, if the 
number of bins is unbounded, the best known asymptotic worst case ratio is 2.1439 [9 . 
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There are also some research on the Multiple Knapsack Problem. Kellerer [T4], first gave a PTAS 
for a special case of this problem in which all the knapsacks have identical capacity [13] • Chekuri and 
Khanna |3J obtained a PTAS for the general multiple knapsack problem. For packing rectangles into 
multiple identical rectangular bins, Fishkin et al. |7J gave a 2 + e approximation algorithm. 

Main results and Techniques: We first observe that the techniques used in Mutilple Knapsack 
Problem [3] are useful for the problem of packing weighted squares into rectangular bins (the bins 
may have different dimensions) with large resources, where large resource means that the height of 
a bin is much larger than the width, and give a Polynomial Time Approximation Scheme (PTAS) 
for the above problem, then apply the PTAS to the problem of packing squares with profits into a 
rectangular bin and get a | + e approximation algorithm. For packing squaures into a rectangular bin, 
we first introduce a simple algorithm by the techniques used in 2D bin packing pQ problems and show 
that if there are too many large squares packed or the rest area for packing squares is not small, the 
algorithm has a nice performance, then we focus on the case in which there are a few large squares 
packed and the rest area for small squares is also few. We propose a novel approach of packing a few 
large items such that the packing does not affect the future small items packing too much, and call 
it corner packing. For packing small squares into the rectilinear polygons which is generated after 
packing large squares into the bin, 

• we first dissect the polygons into rectangular blocks such that the optimal value of packing small 
squares into the blocks is near the optimal value of packing small squares into the polygons, 

• then call the method used in Multiple Knapsack Problem [3], to guess one sublist which has a 
feasible packing and profit at least (1 — e)OPTf J , where OPT^ is the optimal value for packing 
small squares into the blocks, 

• lastly, we exploit the techniques used in strip packing [7J [15] to pack items in each block. 

Worst Case Ratio: We adopt the standard measure worst case ratio to evaluate approximation 
algorithms. For any input list L, let A(L) be the total profit packed by approximation algorithm A 
and OPT(L) be the optimal value. The worst case ratio of algorithm A is thus defined as 

OPT(L) 

RA = s T^(Lr- 

p(-), w(-): Given a square q, we use p(q) and w(q) to denote its profit and area respectively. And 
given a list of squares L = (qi, . . . , q n ), we define p(L) = YH=i P(li) an d w{L) = Y^h=i w { ( li)- 



2 Packing Squares into Rectangular Bins with Large Resources 

INSTANCE: Given an input list S of n squares with profits and a set of rectangular bins B = 
(Bi, -E>2, > • • , B c ) where Bi = (wt, hi) and max{wi, hi} > g 6 * -1 for all i, c, e are constants. 

OBJECTIVE: Maximize the total profit packed in B. 

Based on the ideas from the seminal papers [31 [15], we give a PTAS for the above problem. There 
are three steps in the PTAS. We first guess a subset of squares which can be packed into B and 
whose total profit is near the optimal value through the technique of rounding the input instance into 
0(e _2 lnn) classes. Then for each bin, we guess the number of items packed in that bin from each 
class such that our guess is also near the optimal solution, i.e., we do not lose too much profit. After 
matching items into bins, we use the strip packing algorithm to pack items in each bin. 
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Rounding and guessing: Here, we consider square packing. Since there is an natural order relation 
between any two squares, the techniques used in Multiple Knapsack problem [3] are useful for square 
packing too. We first round the instance into a well structured instance which has 0(e _1 Inn) distinct 
profits, and more items in each profit class have at most (e _1 ) distinct sizes (side length). Then we 
select a subset items which can be packed into the bins and has the profit as least (1 — e) time the 
optimal solution. But, if the items are rectangles, we do not have the above result, since there is not 
an order relation between any two rectangles. 

Lemma 1 Given an above instance I = (B, S) with n items, in polynomial time v = n°^^ e3 \ we can 
obtain instances I\, . . . , I v such that 

• Ij = (B, Sj) for 1 < j < v , where Sj is a sublist of list S. 

• For 1 < j < v, items in Sj have 0(e _1 Inn) distinct profits, and more items in each profit class 
have at most (e _1 ) distinct sizes (side length). 

• There is an index j, 1 < j < v, such that Sj has a feasible packing in B and p{Sj) > (1 — 
0{e))OPT{I). 

Proof. We show how to construct instance Ii,...,I v from / = (B, S) such that one of them satisfies 
the conditions. There are four steps, which are basically from [3]. 

• Guess a value O such that (1 - e)OPT <0< OPT. 

• Round down the profits of items into 0(e~ 1 Inn) classes such that < pj < Pj, where pj is 
the jth item's profit and pj is the one after rounding down, where n is the number of items. 

• Guess a set of sublists based the value O and the rounded profits of items such that one of them 
is feasible to pack into the set of bins (blocks) B and its total profit is at least (1 — 0(e))OPT. 

• Using the techniques in bin packing [5], in each distinct profit class, reduce the number of distinct 
sizes into 0(e~ 1 ) such that we lose the profit at most 0(e)OPT. Hence, totally, each of sublists 
has 0(e _2 lnn) distinct size values and profits. 

Next, we give the details for the above four steps. First, we show how to guess O. Given a sufficiently 
small constant e > 0, let p ma x denote the largest value among item profits. We know the optimal 
solution is bounded by n ■ p m ax ■ So we guess O from the set 

{Pmax(l + e)*|0 < i < 1 + lni +e n < 2c' 1 Inn}. 

(1 + hii +e n < 2e _1 Inn follows from ln(l + e) > e — e 2 /2 > e/2.) Therefore, one of the values in the 
above set is guaranteed to satisfy the desired property for O. 

Given a value O such that max{p max , (1 — e)OPT} < O < OPT, then we show how to massage 
the given instance into a more structured one has few distinct profits. 

1. Discard all items with profits at most eO/n. 

2. Consider the other items and divide all profits by eO/n such that after scaling each profit is at 
most n/e. 

3. Round down the profits of item to the nearest power of (1 + e). 
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It is easily seen that only an 0(e) fraction of the optimal profit lost by the above transformation. 
Since (1 + e) % < n/e, we have 

i < 2e _1 lnn/e < 4e _1 Inn. 

The last inequality follows from n/e < n 2 . Therefore, we can transform the instance into a new 
instance with 0(e _1 Inn) distinct profits such that only an 0(e) fraction of the optimal profit is lost. 

Next we show how to guess the items to pack on the instance with 0(e~ 1 Inn) distinct profits. Let 
h < 4e _1 Inn + 1 be the number of distinct profits in our new instance. We partition the input set 
of squares S into h set Si, Sh with items in each set having the same profit. Let U be the items 
chosen in some optimal solution and let U{ = SiClU. Recall that we have an estimate O of the optimal 
value. If p(Ui) < eO/h, we ignore the set Si; no significant profit is lost. Hence we can assume that 
eO/h < p(Ui) < O and approximately guess the value p(Ui) for 1 < % < h, where P{Ui) is the total 
profit in Ui. More precisely, for each i we guess a value ki £ [h/e 2 ] such that 

h(e 2 0/h) < p(U t ) < fa + l)(e 2 0/h), 

where [h/e 2 ] stands for the set of integers 0, 1,...., [h/e 2 \. 

We show how the numbers k{ enable us to identify the items to pack and then show how the 
values ki,....,kh can be guessed in polynomial time. Given the value ki we order the items in Si in 
increasing size values (side length). Let a% denote the profit of an item in Si. If cij < eO/h, pick 
the largest number of item from this ordered set whose cumulative profit does not exceed ki(e 2 0/h). 
If ctj > eO/h we pick the smallest number of items, again in increasing order of side lengths, whose 
cumulative profits exceeds ki(e 2 /h). The choice of items is thus completely determined by the choice 
of the ki. For a tuple of values k±, kh, let U(ki, kh) denote the set of items packed as described 
above. 

From the above selection, there exists a valid tuple (ki,...,kh) with each ki 6 [/i/e 2 ] such that 
U (ki, kh) has a feasible packing in B and p(U (ki, kh)) > (1 — e)0. 

Now we show that the values ki,....,kh can be guessed in polynomial time. Before that, we 
introduce a useful claim. 

Claim 1 [3] Let f be the number of g-tuples of non-negative integers such that the sum of tuple 
coordinates is equal to d. Then f = ■ If d + g < ag then f = 0(e ag ). 

By Claim[H the number of /i-tuples (k\, kh) with ki 6 [h/e 2 ] and Yli ki — h/e 2 is 0(n°( e 3 )). 

Next we show how to reduce the number of distinct sizes (side length) in each profit class. The 
basic idea is the one used in approximation schemes for bin packing [5j. Let A be a set of g items 
with identical profit. We order items in A in non-decreasing sizes and divide them into t = (1 + 1/e) 
groups A\,...,At with Ai, At—i containing [g/t\ items each and At containing (g modi) items. 
We discard the items in At-i and for i < t — 1 we increase the size of every item in Ai to the size of 
the smallest item in Ai + \. Since A is ordered by size, no item in Ai is larger than the smallest item 
in Ai + \ for each 1 < i < t. It is easy to see that if A has a feasible packing then the modified instance 
also has a feasible packing. We discard at most an e fraction of the profit and the modified sizes have 
at most 2/e distinct values. Applying this to each profit class we obtain an instance with 0(e _2 lnn) 
distinct size values. 

Hence, we have this lemma. □ 

Distributing the selected items into each bin 

After guessing a polynomial number of sublists, next we consider how to distribute the selected 
items in each sublist into bins. Easily to see, the possibilities to assign the selected items into bins is 
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bounded by c n , which is an exponential size of n, where c is the number of bins and n is the number 
of items to be packed. But we can guess a subset from the selected items in a polynomial time such 
that the total profit in the subset is near the optimal solution. 

After step 1, we have (e -2 Inn) classes in the input instance. Let fc, be the number of items of the 
ith class and let if be the number of items of the ith class packed in the jth bins, where 1 < j < m. 

Lemma 2 We can guess a set of numbers h\ in polynomial time such that (1 — e)l\ <h\ < If , where 
1 < e~ 2 In n and 1 < j < c and c is the number of bins. 

Proof. For the jth bin, we guess h\ items from the ith class. If k% < e (i+ £ ) then we can guess a number 

h\ such that h\ = l\ in 0{ ^ + ^ ) time. Else, we guess a number h\ from the set {[(1 + e) x ■ \x = 

1,2, ...} such that {l-e)l\ <h\ <l\. Since h\ < ki, the number of guesses required to obtain a single 
h\ is bounded by g = log 1+e c/e < 0(e -2 lnc), for each class, the total number of guesses for all h\ is 
bounded by g c < 0(e _2c c c ), where 1 < j < c. Therefore for all the 0(e _2 lnn) size classes the total 
number of guesses for is bounded by n e , which is a polynomial of n, where c and e are constants. □ 

Since all the items in each size class have the same profit and by Lemma [2] we have hf > (1 — e)q, 
there exists one assignment which is feasible to B and keeps at least (1 — e) times the profits. Next 
we consider how to packing items into each bin. 

Packing each sublist into each bin: In each bin, we have the following property 

max{w, h} _ 4 
min{u>, h} ~ 

then the techniques used in [151 E! are helpful to pack all squares into the bins. 

First, we first give an important lemma for packing squares into a bin with large resource, called 
cutting technique. 

Lemma 3 Given an input list L of squares with sides at most e and two rectangular bins B\ = (1, a), 
B 2 = (l + 2e,a), then 

(1 - 4e) • OPT(L, B 2 ) < OPT{L, Bi), 
where OPT(L, B) is the optimal value for packing list L into bin B. 

Proof. Now we construct a packing in bin B\ from an optimal packing in bin B 2 and prove its profit 
is at least (1 - 4e) • OPT{L,B 2 ). 

Consider an optimal packing in bin B 2 , we cut B 2 into [^J pieces of slices, say S\, S 2 , . . . , S^j_j 
respectively, such that every slice has an exact width 4e (except the last one), shown as fig. [TJ (Note 
that some squares may be cut into two parts, one part in Si and another part in Sj+i). Then we find 
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Figure 1: Cutting bin B 2 into slices 



a slice Si such that p(Si) < 4eOPT(L, B 2 ) and remove all squares completely contained in slice Si if 
any. Observe that after the above removal, all squares remaining in bin B 2 can be packed into B\. 
Hence, OPT(L, B x ) > (1 - 4e) • OPT{L, B 2 ). □ 
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Lemma 4 1 15\j There is an algorithm A which, given a list L of n square and a positive e, produces 
a packing of L in a strip of width 1 and height A(L) such that A(L) < (1 + e)Opt(L) + 0(l/e 2 ). 

Lemma 5 For packing small squares into a constant number of bins, for each bin, if ™in{iu'fc) — e_4 ' 
then there is a polynomial time algorithm with an output at least (1 — 0(e))OPT , where OPT is the 
optimal value. 

Proof. Given an instance / = (B,S), where S is the set of small squares with profits and B = 
(Pi, B2, ■ ■ ■ , B c ) is the set of rectangular bins, by the method in Lemma [fl we guess a subset Sj C S 
such that Sj has a feasible packing in B and p(Sj) > (1 — 0(e))OPT(I). 

In an instance I p = (B, S p ), we first guess h\ by Lemma [2j Then according to h\ value we assign 
the items to each bin and use the APTAS in Lemma H] to pack items in each bin, where 1 < j < c. 
If in each bin (w, h) the height used by the APTAS in Lemma d] is bounded by (1 + e) max{w, h} + 
0(min{u>, h}/e 2 ), then we keep the assignment otherwise reject the assignment. Since there is a 
S p such that S p has a feasible packing in B. After all the guesses, there is at least one assignment 
remained. For the assignment, we apply the APTAS in Lemma[5]and the cutting techniques in Lemma 
[3] such that in each bin, the profit keeps at least (1 — e) times the optimal value. 

Hence we have this lemma. □ 



3 Previous Algorithms for Packing Squares 

Based on previous techniques used for 2D packing problem [TJ [TT] and the greedy packing (which is 
given in appendix) , we introduce a simple algorithm A\ which is implied in [Tj [TTJ for packing a set 
of squares into a bin (l,h), where h > 1. There are two steps in A\: first group all squares by their 
sizes and guess one group which does not significantly affect the optimal packing and delete it from 
the input list, then pack large items by enumeration, lastly append small items in the "gap" of the 
bin. Next, we give the details of the two steps. 

Grouping: For an integer k = [-], where e < (2h + 2/i 2 ) -1 is sufficiently small and h > 1 is the bin 
height, we select k points in the region (0,1], Pi, . . . , Pk as follows 

Pi = e 6 ' and 1 < i < k. 

Then the interval (0,1] is divided into k + 1 intervals, , ifc+i, where Ii = (Pj,Pj_i], 2 < i < k, 

h = (Pl,l] and/ fc+1 = (0,P fe ]. 

Notation: In the following, given a list L of squares, Lj denotes the list in which all square's sides 
are in interval Jj, w(Li) denotes the total area of Li, p(Lj) the total profits of Li, |Lj| the number of 
squares in Li, where 1 < i < k + 1. 
Packing: 

1. Guess an index i such that OPT(L - L { ) > (1 - e)OPT(L). 

2. Get all feasible packing for Lj_i U • • • U L\, pack each of them into the bin, then partition the 
uncovered space into rectangular bins(blocks) in the method [I] and append Lj+i U • • • U Lk+i 
into these bins by the Greedy algorithm. 

3. Output the one with the largest profit. 

Since there are k + 1 sublists L\, . . . , in L, then the guess in step 1 of A\ is always feasible, 
where k = [-]. After selecting an index i, we define all items in Lj_i U • • • U L\ as large items and 
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all items in Lj+i U • • • U Lk+i as small items. Note that if i = 1 then there are no large items, and if 
i = k + 1 then there are no small items. 

A\s worst case ratio is related to the number of large items in L op t and the rest area for small 
items, where L op t is a sublist of L to produce an optimal solution. 

Fact 1 JJJ/ Given large items with sides larger than e 6 * 1 which can be packed in the bin (l,h), and 
small items with sides at most e 6 ', where i > 1, if the total area of all the squares is at most h — 
£ 4x6 l -l^ y^ ere 2h(l + h) < e _1 , then all can be packed in the bin. 

Lemma 6 J7] [17]/ After packing large items, if the rest area in the bin is at least e~ 1 5, then A\(L) > 
(1 - 2e)OPT(L), where 5 = e 4 * 6 " 1 " 1 . 

Lemma 7 /77]/ Let m be the number of large items in L op t- i) A\(L) > (1 — e)OPT(L) if m = 0; iij 
e fa e A(i)>^(l-e)OPT(I). 

Lemma 8 J7JJ/ Algorithm A± is ran in polynomial time of n. 

4 Corner packing 

To pack squares into a rectangular bin, there are a lot of approaches, the most two studies are NFDH 
PQ and BL. In this section, we first give a new approach, called Corner packing, which includes the 
above two approaches. Then we analyze the corner packing later and show that it is one of the key 
points for improving the worse case ratio. 

During packing squares into the rectangular bin, the uncovered space of the bin may get into the 
rectilinear polygons. The corner packing (shown as in Fig. EJb)) can be regarded as a sequence of 
packing. Every time when one square is packed into the bin, we obey the following rules: 

• select one vertex of the current rectilinear polygons at which the interior angle is 90 degrees, 

• place the square such that one of its corners coincides with the vertex we selected. After packing, 
we get the new rectilinear polygons. 

Note that both NFDH and BL [12] belong to Corner packing, where BL packing is to pack squares in 
the bin as bottom as possible then as left as possible. 

Lemma 9 Assume n squares are packed in the bin by corner packing, then 

i) there are at most 4 + 2n vertices of all the rectilinear polygons, 

ii) there are at most 2 n (n + 1)! possibilities to pack these n squares in the bin by corner packing. 
(refer to the proof in the appendix). 

5 A Refined Algorithm for Packing into a Rectangular Bin 

Let m be the number of large items in L opt , where L opt is a sublist of L to get the optimal solution. 
By Lemma [7] if m is very large, then algorithm A± has a good performance. So, to improve algorithm 
Ai, we have to study the case in which m is very small. Note that when the bin is a unit square, the 
situation becomes a relatively simple. Since when m = 1 we can transform the original packing into 
a special strip packing; when m = 2, 3 we can estimate there must be much more space left for small 
squares than the wasted area. This is the main idea in Harren's paper [TT]. If the bin is not longer a 
unit square, his algorithm does not work very well. To improve algorithm Ai, we are faced with two 
problems: 
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• How to pack a few large items such that the packing does not affect too much the future small 
items packing? (how to allocate large items in the bin.) 

• How to pack small items in the gaps (rectilinear polygons) generated after packing large items? 

Next, we give our solutions for the above questions and propose a refined algorithm called A% with 
the worst case ratio | + e. 

5.1 Packing a few large items 

Recall that, given an i > 1, if a square's side length is at least e 6 ' 1 then it is called large, else its side 
length is at most e 6 ' then it is called small. And there is a gap between large items and small items, 
which is very important for packing large items. 

Next we show that corner packing is a good packing which does not significantly affect future small 
square packing when there are a few large items. 

Lemma 10 Let m be the number of large items in L op t. If m < 4, then (1 — e)OPT(L) < OPT(L, *), 
where OPL(L) is the optimal value of packing L into the bin and OPT(L, *) is the optimal value of 
packing L into the bin such that all large squares are packed by corner packing. 

Proof. Here, we just give the proof when m = 4, since the proof for m = 1, 2, 3 is involved. Let a, b, c, d 
be the four large squares in an optimal packing L opt . Without loss of generality assume a,b,c,d are 
placed in the bin as Fig. [2](1) . Note that a large item has side at least e 6 * and a small item has 
side at most e 6 ' , where i > 1. We cut the bin into three parts, two rectangular blocks / = (xi,yi), 
= (x2,U2) and a rectilinear polygon P as shown as Fig. [2]J2). Now we define two new rectangular 
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(1) (2) ( 3 ) 

Figure 2: An optimal packing vs. its corner packing 



blocks I e = (xi + 2e 6 \yi) and II t = (x 2 ,y 2 + 2e 6 '). Then we have 

OPT{L,B) <OPT{L,I e UlI e UP), (1) 

where OPT(L, B) is the optimal value of packing L into the bin B and OPT(L, I e U -ZT e U P) is the 
optimal value of packing L into three rectilinear polygons I e U II e U P. This can be seen as follows, 
all squares packed into the bin B as shown in Fig. [2|L) can be packed into three rectilinear polygons 
I t U II e U P. By Lemma [3l for any list L of small squares, we have 

(1 -Ae 61 )OPT(L,I e ) < OPT(L,I) and (1 - 4e & )OPT(L, II e ) < OPT(L,II). 

Then 

(1 - 4e & )OPT(L,I e U II e U P) < 0PT(L,IUIIUP). (2) 
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So, by ©, ©, 

(1 - Ae & )OPT(L, B) < OPT(L, IUIIU P). 

And we have OPT(L, I U II UP) < OPT(L, *a), where OPT(L, *a) is the optimal value of packing 
L into the bin with a at one corner, shown as in Fig. [2^3). Hence 

(1 - 4e e ')OPT(L, B) < OPT(L, *a) . (3) 

By the similar proof, we have 

(1 - 4e 6l )OPT(L, *o) < OPT(L, *ab) 

and 

(1 - Ae & )OPT{L, *ab) < OPT{L, *abc), (1 - Ae & )OPT{L, *abc) < OPT{L, *abcd) 

where OPT(L, *X) is the optimal value of packing a set X into the bin with all the items in X at 
corners of the bin. Therefore, we have 

OPT(L, *abcd) > (1 - 4e 6l ) 4 OPT(L, B) > (1 - e)OPT{L, B). 

The last inequality follows from e < 1/2. 

Hence, this lemma holds. □ 



5.2 Packing small items into rectilinear polygons 

After packing large items, the uncovered space in the bin may be a set of rectilinear polygons. Our 
strategy for packing small squares into the polygons are below: 

• Dissect the polygons into rectangular blocks such that the optimal value of packing small squares 
into the blocks is at least (l — e)OPT p , where OPT p is the optimal value for packing small squares 
into the polygons. 

• To pack small items into blocks, we use the PTAS in Section [2] of packing squares into rectangular 
bins with large resources. 

Dissection: After packing few large squares into the bin by corner packing, we dissect the recti- 
linear polygons into rectangular blocks, such that the dissection does not affect the optimal packing 
insignificantly. 

Lemma 11 // there are at most 4 large squares packed, and the total area of the large squares packed 
is at least h — e 4xgl ~ 2 , then there exist a dissection of the polygons (the uncovered space of the bin) 
into blocks such that 

OPT b > (1 - e)OPT p , 

where OPTf, (OPT p ) is the optimal value of packing small squares into blocks (polygons). 

Proof. In this proof, we just give our dissection for four large squares packed, shown as in Fig. [3] and 
H] (by dotted lines), since the number of large squares is less than 4, we have the similar dissection. 

Except for case (d), if the rectilinear polygons are dissected into blocks as shown in Fig. [21 then 
we have a set of five rectangular blocks B = {Bi}, where Bi = (wi,hi) and W{ < hi for 1 < i < 5. 
(Otherwise we can exchange Wi and hi) Using the same techniques in Lemma [T0| we define a new set 
of five blocks B = {-E^}, where B i = (u>i,hi + 2e 6 ') for 1 < i < 5. Let P be the polygon(s) after 
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(a) (b) (c) (d) (e) 

Figure 3: Possible packing 
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S4 
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Si 


Bi 




S2 



S3 



B3 S4 



B2 : ^> 



Si 



B4 



Bi S2 



(dl) (d2) 

Figure 4: Two possibilities for the case (d) 



packing large squares in the bin, and L be a list of small squares. Since each small square has side 
length at most e 6 \ we have 

OPT(L,B') > OPT(L,P) > OPT(L,B), 

where OPT(L,X) is the optimal value for packing L into X. Since each large square has side length 
at least e 6 ' and hi > Wi, we have hi > e 6 * . By Lemma we have 

OPT(L,B) > (1 - e)OPT(L,B'). 

Hence we have OPT b = OPT(L, B) > (1 - e)OPT(L, P) = (1 - e)OPT p . 

Next, we study the case (d) of Fig.[3]and prove that our strategy shown in Fig. Estill works. There 
are two possibilities for the case (d). We assign Si, S2, S3, S4 to the four large squares as shown in 
Fig. [U where Si = (sj, Si). And the polygon is dissected into 5 blocks B\, . . . , B§, where Bi = (wi, hi). 

From our dissections in Fig. [H (by dotted lines), we have 

max{u>j, hi] > e 6 , for 1 < i < 4. 

And since the total area of the large items in the bin is at least h — e 4x61 1_2 ) i.e., the total area of 
the blocks is at most e 4x6 * ~ 2 . We have 

£ 4x6 l ~ 1 -2 

min{wi,hi} < — — < e 2x6 ' . 

e° 

The last inequality follows from i>2. 
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Let 5 = e & , next we prove that < 5 2 in the cases (dl) and (d2). It is trivial to see w$ < 5 2 
in the case (dl), since < w\ = min-fiui, hi}. Now, we consider the case (d2), since s\ + S2 < 1 and 
si — S2 < I14 < 5 2 , we have 

1 + 5 2 
81 ■ 

Since S3 + S4 < 1 and 54 — S3 < hi < 5 2 , we have 

1 + 5 2 

So, 

w 5 = si + s 4 - 1 < 1 + S 2 - 1 = 5 2 . 

Therefore, we have in the cases (dl) and (d2), W5 < 5 2 < ew^ and h§ < /14, i.e., to compare with 
block I?4, block is sufficiently small and can be ignored. So 

OPT(S, B~) > (1 - 0(e))OPT(S, B), 

where S is the set of small squares, B~ = \jf =1 Bi, B = U^ =1 Bi. Hence when we pack small squares 
into blocks Ll| =1 -Bj, we just consider L)f =1 Bi. Then by the similar proof for other cases, we have 
OPT b > (1 - 0(e))OPT p . 

Hence, this lemma holds. □ 



5.3 Algorithm A2 and its analysis 

Next, we give the details of algorithm A2 which is similar to A±. 
Description of Algorithm A2 

1. Group items and guess an index i such that OPT(L — Li) > (1 — e)OPT(L) and divide the 
remaining into two classes, say large and small, 

2. Get all feasible packing of Lj_i U • • • U Li, for each of them, 

(a) if there are at least 4 items or the total area of items is at most h — e 4x61 1 ~ 2 , then pack 
large and small squares by algorithm A\. 

(b) else locate large items as Fig. [3] and divide the gaps into blocks as Fig. [3] and [H lastly apply 
the method in Lemma [5] for small items. 

3. Output the one with the largest profit. 

Theorem 1 For any input list L, °^^ > < |(1 + 0(e)), where e is sufficiently small. 

Proof. To consider an optimal packing solution L op t, if there are at least 5 large items in L op t or the 
total area of large items in L op t at most h — e 4x61 1_2 5 by Fact[lJ Lemmas [6] and [71 

A 2 (L)>l(l-2e)OPT(L). 

Else, the total area of the large items in the bin is at least h — g 4x6 ' 1 ~ 2 and there are at most 4 
packed. 
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By the dissection of the polygons into rectangular blocks, shown as in Fig. [3]and[U in each block 
(w, h), we make sure that max{ro, h} > e 6 * . So, 



max{w, h} maxjw, h} £ 2x6 l 1 ^ 
mm{w.h} ~ £ 4x6'-i-2~ - F 4x6'- 1 -2 - e 

The last inequality follows from i > 2. (remember when i = 1, there is no large item.) By Lemmas 
[TUl 03 we have 

A 2 (L) > (1 - 0{e))OPT{L) 

By Lemma the time complexity of Algorithm A2 is a polynomial time of n. Hence, this theorem 
holds. □ 

6 Concluding remarks 

Note that algorithm A2 can be expended to multi-dimensional cube packing. 
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7 Appendix 



7.1 The proof for Lemma [9] 

Proof. We use induction to prove part i). When n = 0, there are 4 vertices in the bin. When n = 1, 
there are at most 6 vertices in the rectilinear polygon. So, we assume that when n = k, part i) holds, 
i.e., after packing k items in the bin, there are at most 4 + 2k vertices in the rectilinear polygons. 
When n = k + 1, we use one of 4 + 2k vertices and generate at most 3 vertices, hence the total number 
of vertices is at most 

4 + 2/c-l + 3 = 4 + 2(A; + l). 

Then we can see there are at most 4 + 2(i — 1) ways to pack the i-th square, where i > 1. Hence to 
pack n items in the bin, there are at most 

n n 

[](4 + 2(i - 1)) = 2 n \\(i + 1) = 2 n (n + 1)! 
i=l i=i 

possibilities. □ 



7.2 NFDH packing 

NFDH (Next Fit Decreasing Height) [18J. NFDH packing behaves as follows: First sort all squares by 
their heights, then pack them in the bin from the largest one level by level as shown in Fig. [5ja). In 
each level, pack them by Next Fit, namely, if the current level cannot accommodate the next item, 
then open a new with height equal to the current item's height. We repeat this procedure, until there 
is no space for a new level in the bin. Here is a key property of NFDH. 



level 3 



level 2 



level 1 



(a) (b) 




Figure 5: NFDH and Corner packing 



Greedy Algorithm 

1. Sort the input list L such that 444 > > 444- 

r w(Ai) — — w(A k ) 

2. For i from 1 to n do 

if (o,i > e) and (pi > e) then 

(a) Find a maximal index m such that (A\,A2, ...,A m ) can be packed into the current bin by 
NFDH and pack (A\, A2, A m ) into (cij, b{). 

(b) Then update list L and re-index L. If L becomes empty then finish packing. 
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